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Abstract: In this paper we analyze an exactly solvable model consisting of an inertial 
Unruh-DeWitt detector which interacts linearly with a massless quantum field in Minkowski 
spacetime with a perfectly reflecting flat plane boundary. Firstly a set of coupled equations 
for the detector's and the field's Heisenberg operators are derived. Then we introduce the 
linear entropy as a measure of entanglement between the detector and the quantum field 
under mirror reflection, and solve the early-time detector-field entanglement dynamics. Af- 
ter coarse-graining the field, the dynamics of the detector's internal degree of freedom is 
described by a quantum Langevin equation, where the dissipation and noise kernels respec- 
tively correspond to the retarded Green's functions and Hadamard elementary functions 
of the free quantum field in a half space. At late times when the combined system is in a 
stationary state, we obtain exact expressions for the detector's covariance matrix and show 
that the detector-field entanglement decreases for smaller separation between the detector 
and the mirror. We explain the behavior of detector-field entanglement qualitatively with 
the help of a detector's mirror image, compare them with the case of two real detectors 
and explain the differences. 
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1 Introduction 

The effects on a quantum field due to the presence of boundaries (such as a mirror or 
a dielectric plate) has been studied a long time ago, as in the famous Casimir effect [1]. 
Moving mirrors in a quantum field have been used as analog models [2] of Hawking and 
Unruh effects [3, 4]. Quantum field theory in spacetimes with nontrivial topology [5] has 
also been studied since the seventies. Today these subject matters fall under the broadened 
scope known as quantum field theory under external conditions [6]. 

Quantum entanglement is, according to Schrodinger, "t/ze characteristic trait of quan- 
tum mechanics, the one that enforces its entire departure from classical line of thought" [7]. 
The nature and behavior, even the definition for more than two parties, of entanglement 
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has attracted much attention in the last two decades for both theoretical and practical rea- 
sons, the former because it is the cornerstone issue at the foundation of quantum mechanics 
and the latter connected with the rapid development of quantum information sciences. 

Quantum entanglement of systems composed of detectors ^ and mirrors, coexisting 
with and/or mediated by a quantum field is also receiving increased attention as an ac- 
tive research branch in a newly emergent field known as relativistic quantum information 
[9] - relativistic because the effects of the quantum field is highlighted where traditional 
treatments based on nonrelativistic quantum mechanics prove inadequate. (For a recent 
review of problems in RQI involving UD detectors see, e.g., [10]). An important theoretical 
issue which quantum field theory brings forth in contrast to quantum mechanics is quan- 
tum nonlocality, a vivid illustrative example of this is given in [11], where entanglement 
between two causally disconnected objects can be created by their local couplings with a 
common quantum field. 

Amongst theoretical problems on quantum entanglement of current experimental in- 
terest [12] we mention two aspects: 1) Time evolution, or entanglement dynamics: Even 
in the simplest system of two two-level (2LA) atoms interacting via a common quantum 
field already one sees a diverse spectrum of interesting behavior, ranging from sudden 
death [13], touch of death, revival, to staying always alive, and features such as dynami- 
cal generation, protection, and transfer of entanglement between subsystems [15]. Model 
studies of oscillator-field entanglement [14] with experimental ventures have been carried 
out [16] as well as mirror-field entanglement in conjunction with measurements in LIGO 
detectors [17], the latter also serving as preparatory studies for quantum superpositions of 
mirrors [18] in macroscopic quantum phenomena. 2) Spatial dependence: Entanglement 
not only changes in time but also depends on spatial separation, as shown in model studies 
for 2LAs [19-21] and harmonic oscillators [22] interacting with a common quantum field. 
Practical application of these studies span from the design of quantum gates to quantum 
teleportation [23]. 

In this paper we study the model consisting of an inertial UD detector which interacts 
linearly with a massless quantum field, the latter being restricted to a half space by a per- 
fectly refiecting boundary, or equivalently, a perfect mirror. The presence of the mirrors 
is known to influence radiative properties of physical systems, such as the rate of sponta- 
neous emission of excited atoms [24]. Here we probe such systems at a deeper level, asking 
how the mirror affects the delicate property of entanglement between the detector and the 
quantum field at zero temperature. 

The effect of the mirror can be understood from two perspectives. Prom a purely 
technical consideration, the presence of the mirror breaks Poincare symmetry and thus 
modifies the spectral density of the field, which accordingly alters the entanglement and 
its dynamics of the detector-field system. From physical considerations one can exploit the 
symmetry in the set-up and see if this problem of a detector interacting with a mirror- 
modified field can be solved more simply with the use of a mirror image, as is familiar in 

^By detector we refer to a physical entity with some internal degrees of freedom, e.g. an atom, an 
electron with spin, or an oscillator. We will be working with the Unruh-DeWitt (UD) detector specifically 
[4, 8]. 
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classical electromagnetics. From this perspective, one may even ask whether a detector 
can be entangled with its own mirror image, leading to the notion of self-entanglement ( 
autangle or ipso-tangle ). We will pursue this problem in a cautious manner, i.e., carry out 
the full calculation without assuming any mirror image; only after we obtain the results 
will we then ask if one could interpret them in a mirror image way and answer the query, 
"Can one get entangled with one's mirror image?" 

The paper is organized as follows: In Section 2 we set up the problem and derive a 
formal set of equations for the Heisenberg operators of the detector's internal degrees of 
freedom and the quantum field. We impose the boundary conditions introduced by the 
presence of the mirror, then calculate the retarded Green's functions and Hadamard ele- 
mentary functions of the altered field configurations. In Section 3 we introduce the linear 
entropy as a measure of quantum entanglement and calculate the early-time dynamics of 
entanglement between the detector and the field in our setup. We show how the entangle- 
ment between the detector and field in a half space under mirror reflection evolves in time, 
and how it depends on the distance between the detector and the mirror. In Section 4 we 
derive a quantum Langevin equation for the detector's internal degree of freedom includ- 
ing the back-reaction of the altered field. We seek late-time solutions to this equation and 
calculate the covariance matrix of the detector's canonical variables at late times when the 
combined system is in a stationary state. 

Details are contained in Appendix A. Finally in Section 5 we explain in what sense can 
one describe this situation in terms of entanglement of the detector with its mirror image 
and why it decreases as the detector moves closer to the mirror, a somewhat counterintuitive 
finding. We also compare this situation with the case of two inertial detectors in free space 
(calculations placed in Appendix B) and explain their physical differences. 

2 Unruh-DeWitt Detector in Half Space under Mirror Boundary Con- 
dition 

The gist of the matter for this paper is to quantify the change of entanglement between the 
detector and the field because of the presence of a mirror. To achieve this we first derive 
the Heisenberg equations of motion for both the detector's internal degree of freedom and 
the field under the boundary condition of the mirror. After coarse graining the field, we 
obtain the reduced dynamics of the detector which already contains the back-reaction of 
the field on the detector. Then by solving the reduced dynamics of the detector one can 
obtain its time-dependent covariance matrix. This covariance matrix describes the mixed 
state of the detector's internal degree of freedom, according to which the detector-field 
entanglement can be quantified because the combined detector-field system is in a pure 
state. 

2.1 Dynamics of Detector-Field System with Mirror 

Consider an UD detector located at position at a vertical distance of L/2 from the 
mirror plane at X3 = 0. The action of the total system, which describes the detector with 
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internal degree of freedom Q interacting linearly with a massless scalar field <I> in a half 
space defined by X3 > with coupling constant Xq, is given by 

S =^Mg j dt {Q"^ -nlQ^) + ^ j dt d^xdf,^d^^ + Xq j dtQ{t)^{-iLq,t), (2.1) 

where Q acts like a harmonic oscillator with mass Mq and bare natural frequency f^o- Due 
to the linearity of our system the equations of motion for Heisenberg operators (carrying 
hats), 

MqQ{t) + Mq^fQ{t) = XqH-^q, t) , (2.2) 

□|.(x,t) = Ag<53(x-Xg)Q(t), (2.3) 
i(x||,X3 = 0,0 = 0, (2.4) 

have the same form as the equations for the corresponding classical variables. Equation 
(2.4) specifies the Dirichlet boundary conditions imposed on the field where the mirror 
surface is located, namely, in the X3 = plane and xy = (xi,X2). To obtain the back- 
reaction of the field on the detector we first solve (2.3) and then plug the solution into 
(2.2). The solution for ^>(x, t) is given by 

6(x, t) = $o(x, t) + A, /* dt'G%,it, x; t', ^q)Q{t') (2.5) 

Jt, 

where $o(x, t) is the homogeneous solution to equation (2.3) which describes the dynamics 
of the source-free field (without Q) in the presence of the mirror and G*g^(x, y, t') is the 
retarded Green's function for the field. Plugging the solution for the field operator into 
the equation of motion for the detector we obtain the following equation governing the 
dynamics of the detector with the effects of the field already incorporated, 

Mq^{t) + Mq^'Qit) - Xl f dt'G%t{t, X,; t' , x,)Q(t') = A,io(x,, t). (2.6) 

Jti 

This is what we mean by the 'field-influenced' dynamics of the detector. The term con- 
taining the retarded Green's function describes how the detector transfers energy to the 
field, as will be explained in further detail later, and the right hand side acts similarly to 
a Langevin forcing term describing how quantum field fluctuations drive the detector. 

2.2 Green's Functions of Free Field in Half-space 

In order to solve the equation of motion for the detector's internal degree of freedom (2.6), 
we need to know the evolution of free field in a half space. The free field operator $0 which 
satisfies the Klein-Gordon equation and vanishes in the X3 = plane can be decomposed 
into normal modes: 

l>o(x,t)=/" d^k^-i^e-''^^+'^r^\\smk3X3b]^ + H.c. (2.7) 

Jk3>0 V ^TT'^UJ 

where cj = |k|. Notice that the normalization factor is different compared to the case of 
free space. 
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Writing no(x', t) = '&o(x', t) and demanding that the equal-time commutation relations 
are satisfied in half-space, that is, 



[clo(x,t),no(x',t)] =i<5(3)(x-x'), (2.8) 

we can infer that the mode expansion coefficients can be identified with the creation and 
annihilation operators 

[^k,^L] =^^1^1^', [k,k]=0, [bl,bi] = 0. (2.9) 

In the above commutators the vectors are restricted to a half space, i.e. X3 > and > 0. 

With the expression for the free field operator in hand we now calculate the retarded 
Green's function and Hadamard elementary function. The retarded Green's function quan- 
tifies the field sourced by the detector and describes how energy is dissipated from the 
detector into the field: 

G*,(x,t;y,0 ^ i9{t-t')[^o{^,t),^o{y,t')] 

= e{t-t')^ [ d^k-sm{u{t-t'))e^^\\<''\\-y\\hm{k3X3)smik3y3) 

/•oo 

= B{t-t')j dojsm{u{t-t')) ■ I{oj;yL,y), (2.10) 
Jo 

where the spectral density I{uj; x, y) is given by 

/(a;;x,y) = ^/ dO e*''l|-(''ll~yil) sin(w cos 6'x3) sin(u; cos ^a). (2.11) 

27r'^ Jq Jo 

The Hadamard elementary function, on the other hand, describes the quantum fluctuations 
of the field. It is given formally by the anti-commutator of the field operator: 



^ / d^k-cos{uj{t-t'))e^^\\<''\\-y\0sm{k3X3)sm{k3y3) 



G%{^,t;y,t') ^ ({l>o(x,t),$o(y,0}) 
_ 1 

POO 

= / dujcos{uj{t-t')) ■ I{uj;yL,y), (2.12) 
Jo 

which is derived from the same spectral density (2.11). 



3 Early Time Dynamics of Detector-Field Entanglement 

Understanding how entanglement in various types of quantum systems evolves in time 
is of basic importance for the purpose of quantum information processing [12] and can 
provide new insight into certain foundational issues involving apparent information loss and 
unitarity violation. For example, the dynamics of entanglement between an accelerating 
UD detector and a massless field in 3+1 dimensional Minkowski spacetime has been shown 
to have implications for the issue of black hole information [25]. 
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In [22] the dynamics of the entanglement between two inertial detectors with identi- 
cal couplings to a common massless scalar field was investigated and it exhibited several 
interesting features. There are roughly three stages of evolution if the two detectors are 
properly separated and weakly coupled to the field. The first stage is at early times, up 
to t ~ 0(1/Ag) when accumulated effects of field- mediated mutual influences between the 
detectors, which manifest themselves in the field's retarded propogators, are still weak and 
the reduced dynamics of the detectors are dominated by the infiuences of vacuum fluctua- 
tions of the fleld. As the effects of fleld-mediated mutual influences between the detectors 
gradually gain strength, the system will enter an intermediate stage in which the detectors' 
reduced dynamics become quite complicated. In the late-time limit, or the final stage, the 
detector approaches a stationary state in time. 

In particular, in [22] it was shown that at early times, after causal contact has been 
established between two detectors, an oscillatory pattern of entanglement emerges which 
varies with the distance between the detectors. This is not due to mutual infiuences since 
in the weak coupling limit the mutual influence effect is always small compared to effect of 
vacuum fluctuations. 

For our model, we expect a similar partition of the history in three stages. Here in this 
section we study the early-time behavior of our model, assuming that initially the detector 
is in its ground state uncorrelated with the field, which is also in its ground state under an 
external Dirichlet boundary condition. We will show that the detector-field entanglement 
develops a spatial oscillatory pattern characterized by the renormalized frequency of the 
detector fi^ (we adopt the convention c = 1). Moreover, the growth rate of the detector- 
field entanglement also oscillates as a function of the detector-boundary spacing. 

3.1 Measure of Detector-Field Entanglement 

In our setup, the detector and the field together form a closed quantum system which 
undergoes unitary evolution. Since the initial state of the entire system is Gaussian and 
the action is quadratic, the reduced quantum state of the detector at any time will remain 
Gaussian. For such systems the amount of bipartite entanglement between the detector and 
the field can be quantified by the linear entropy of the detector's reduced density matrix 
[26], which takes on values larger than when the detector and field are entangled. Here 
the linear entropy is defined as 5^ = 1 — "P, where V = Trp^ is the purity of the detector's 
reduced quantum state. For our general mixed Gaussian state, the purity is given by 
V = l/(2-v/det V) where V is the detector's covariance matrix defined as Vjj = {Oi,Oj), 
in which (Oj, Oj) = Tr(/5a • {Oi, Oj})/2, and O = {Q, P). Mathematically we know that V 
is always less than or equal to 1. Smaller purity means the reduced density matrix of the 
detector is less pure, and correspondingly the detector is more entangled with the field. 

3.2 Mode Decomposition and Exact Solutions for Mode Functions 

We perform the following mode decompositions for Heisenberg operators of the detector 
Q{t) and the field with respect to the creation and annihilation operators of the 
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detector d, a) and those of the field 6, 

,{t)a + ql{t)a) 



Q{t) 



+ 



h 



/a(x)a + /*(x)a 



+ 



fc3>o V2^V 2tj L 



(7+(t,k)6k + g-(t,k)&T , (3.1 




/+(x;k)5k + /-(x;k)6j^ ,(3.2) 



'fc3>o v27r3 V 2uj 

where is the renormalized internal frequency of the detector. From here on subscript 
'a' is always associated with the detector's internal degree of freedom. 

By plugging 3.1 and 3.2 into 2.2 and 2.3, respectively, we find the following equations 
of motion for the mode functions 



L 

{df + 2^^dt + nl)q+{t,'k) = -'^9{t-L)q+{t-L,k) + ^fo+{t,^,;k), (3.4) 



(3.3) 



in which /o+(t,Xg;k) = e «'^*+««^||-x|| ^nd 7^ = A^/(87rMg). Since we have assumed 

iut „;„ kgL 



iq — ->5 

the detector's position to be x^(f) = = (0, 0, L/2), we have /o+(Xq, t; k) = e"*'^* sin ^ 



According to our initial condition, the solutions have to satisfy the initial conditions 
/+(0,x;k) = e^'^ll-^ll sinA:3X3, 9t/+(0,x;k) = -icoe^r^'W sinfes^a, qM = 1, dtqM = -i^r, 
and /a(0,x) = 9Ja(0,x) = g+(0;k) = dtq+{0;k) = 0. 

The solution to (3.3) can be written as the expansion 



'?a(i) = 5^gi")(t), 



(3.5) 



n=0 



with the n-th order retarded influences gi"^ by the (n — l)-th order backreaction sourced 
from the detector then reflected by the mirror: 



qi^Ht) 



7W, 



(3.6) 



q^\t) = / dnGr{t,n) 



27, 
L 



{t-L)x 
27, 



... J dTnGr{Tn-i " r„) (-^) " nL)#)(r„ - uL) . (3.7) 

The expansion of qa{t) is truncated at n = with the minimal N satisfying NL > t, 
because the least time it takes for the n-th order field-mediated influence to affect the 
detector is nL. Here Gr{t,T) is the retarded Green's function which satisfies {d^ + 2'^qdt + 
nl)Gr{t,T) = 6{t-T) and 

2 ■ ' 



it) 



n 



n 



(3.^ 



is the homogeneous solution satisfying the initial conditions. 

The solution to (3.4) can be written in a similar fashion as above, also truncated at 
nL > t, 

g+(t,k) = J^gf (t,k), (3.9) 



n=0 
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where 



q+\t,k) = J driGr(t,Ti)Ag/o+(Ti,Xg;k) 



sin^ [(Ml - M2)e-*'^* + 



(Mae^™ - Mie"^^*) e'^^*] , (3.10) 



i"-'^(t,k) = y driG,(t,Ti) (^-^^^(ri-L) ^ dT2Gr{n,T2) (-^)^(r2-2L)..-x 

y" dTnGr{Tn-l - Tn)0{Tn - (n - l)L)^/o+ (t„ - (n - l)L,Xq;k) , (3.11) 

with Ml = l/[2{-uj - i-fg + n)] and M2 = l/[2(-a; - ijg - Q)]. 
3.3 Zeroth-Order Correlation Functions 

For the factorized initial state assumed in this paper with both the detector and the field 
being in their ground states and uncorrelated, the elements of the covariance matrix can 
be decomposed into two parts corresponding to the two sets of operators in the mode 
decomposition, namely, 



{Q{t),Q{t)) 



1 

2a 



ka(t)P + 



fc3>0 



d^k 1 

2tt^2uj 



= m),Q{t)), + {Q{t),Qit)), 



(3.12) 



Here 'a' corresponds to the detector and 'v' corresponds to the influence of the field's 
vacuum fluctuations. {P(t), P{t)) and {Q{t), P{t)) are similar. 

In the weak coupling limit, the effect of reflected influences correspond to terms in 
(3.5) and (3.9) which are of higher (than zeroth) order; consequently, at early time (up to 
t 1/75) accumulated effect of reflected influences is always small. Thus for the purpose 
of studying the early time behavior of entanglement, we now ignore the contribution of 
reflected influences and restrict ourselves to the lowest order correlations, which read 



(Q(t),Q(t))f = 



1 



n 



(Q(i),Q(t))f = 



'''' .f(.,k)^ 



fc3>o 27r3 2a; 
d'^k 1 



(3.13) 
(3.14) 



fc3>0 
1 



27r3 2u V'''' ^ ' 1^^"^' " ^2)e-*"* + (Mse*^" - Mie' 



w / sincjL 
M„0; 47r2 / 2 I 



I (Ml - M2)e-^^* + (M2e^^* - Mie~'^^)e-^^H^ 



and so on. 

Notice that the a-parts of the zeroth order correlators (.•.)a'^'' do not depend on the 
distance L/2 between the detector and the boundary, but the v-parts induced by vacuum 
fluctuations do have such dependence. 
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In [22], for the case of two inertial detectors distance L apart, which are located in 
free space with identical couplings to the field, the zeroth order correlators due to vacuum 
fluctuations are given as, for example {Mq = 1 and d in [22] is replaced by L), 

I (Ml - M2)e~*^* + (Mse^^* - Mie"^™) e-T«*f . (3.15) 

Physically, the v-parts of the zeroth order correlators (.••)v''^ effectively measure the 
response of the detector to vacuum fluctuations of the field. The similarity between the 
integrands of {Q{t),Q{t))v in Eq.(3.14) and {QAit),QB{t))v in Eq.(3.15) is not surprising. 
In Appendix B we show that, if we have two inertial detectors C and D in free space 
at a distance L apart, with coupling constants of the same magnitude but in opposite 
signs, then {Qc{i) + QD{t))/2 obeys the same equation of motion as the one for Q(t) in 
our model. Thus we see that the self correlator {Q{t),Q{t))^ here has the same value 
^ HQcQc) + {Qd,Qd))/^ + {QctQd) /"^ at t and contains the part of correlations of 
vacuum fluctuations in free space which is odd with respect to the = plane, whereas 
{Qc{t)-,QD{t))v oc {—Xq)\ here has exactly the same value of — {Q A{t) ■, Q B{t)) v in [22] 
because there the two detectors are identically coupled to the field and so {Qa,Qb)v fx 

XqXq. 

3.4 Early-time Dynamics of Detector-field Entanglement 

With the previous results we can now demonstrate the evolution of detector-field entan- 
glement as the distance between the detector and the boundary changes. Here we study 
the dependence of the linear entropy on L and t numerically, as shown in Figure 1, and 
provide physical explanation for its behavior. 

Our numerical results reveal the following behaviors: 

• For a given L, the entanglement between the detector and the field increases mono- 
tonically at early times, showing that the detector is getting more and more entangled 
with the field after the interaction is turned on. (See Figure 1 (upper row) and (lower 
right).) 

• Within the light cone, at every given time in this stage, the detector-field entangle- 
ment exhibits an oscillatory behavior with spatial frequency being 0,. This shows 
that the growth rate oi Sl with which the detector gets entangled with the field also 
oscillates as a function of L, as shown in Figure 1 (upper row) and (lower left). In [27] 
the transition rate of a detector close to an infinite plane boundary with Neumann 
boundary condition was shown to exhibit similar oscillations. 

Such oscillatory behavior in L is solely due to the oscillatory behavior of the zeroth 
order correlators corresponding to {Rc, Rd)v, R = Q, P, given in the last paragraph of the 
previous subsection. After cos ^ in ^3 = w cos 6 has been integrated over the interval [0, 1] 
in the integration corresponding to {Qc, Qd)v in (3.15), the field modes in resonance with 
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Sl 




Figure 1: Entanglement dynamics at early time where only zeroth order correlations 
contribute. Here 7^ = 0.02 and 17 = 5. The upper left plot shows the linear entropy as a 
function of L and t, with the upper-right plot being its contour plot. The lower- left plot 
shows the dependence of linear entropy on L at a given instant of time whereas the lower 
right plot exhibits how the linear entropy evolves with time for the detector located at a 
certain distance. 



the detector have constructive (destructive) interference at the local maxima (minima) of 
Sl against L at fixed t at early times. 

In [22] the early time dynamics of entanglement between the two detectors exhibit 
similar oscillatory dependence on the distance between them, which is also due to distance- 
dependent correlations of vacuum fluctuations experienced by the detectors. More precisely, 
from (3.15) one can see that the field modes with u = nir/L, n = 1,2,3,---, has no 
contribution at all to the integration of {Qa, Qb)v, and those satisfying tanwL = uL give 
the maximum and minimum values of the factor sina;L/(ti;L) in the integrand. In the 
weak coupling limit the integration in (3.15) is mainly contributed by the poles at a; ~ ±17, 
namely, those mode on resonance with the detectors. So {Qc,Qd)v ~ for L mr/Q, 
while {Qc^Qd)v has local maximum or minimum values when L is about the solution of 
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tanfiL = QL in the weak coupling limit. 



4 Late-Time Stationary Limit of Detector-Field Entanglement 

In the late-time limit, as information flows from the detector to the field and gets dissipated 
through the interaction with the field, the entire system approaches an asymptotically 
stationary state. Under sufficiently general conditions, one can say that, supposing the 
field is initially in the vacuum, the late-time asymptotic detector-field entanglement is 
independent of the initial condition of the detector as long as the detector-field coupling is 
linear. As our result will show, the detector-field entanglement is a delicate manifestation 
of quantum correlations between the detector and the field in the presence of the mirror. 

4.1 Quantum Langevin Equation and Covariance Matrix at Late-times 

One can compare the dynamics of the detector under the infiuence of the field to the well- 
studied quantum Brownian motion (QBM) model, namely, the retarded Green's functions 
and Hadamard elementary functions correspond to the dissipation and noise kernels re- 
spectively and the function I{uj;x,y) as in (2.11) represents the spectral density (see, e.g., 
[32, 33] and the references therein.) In the same vein Eq. (2.6) can equivalently be written 
as a quantum Langevin equation 

M,Q{t) + MqQ^Qit) - f dMt, r) • Q(r) = i{t), (4.1) 

•Jti 

in which the dissipation kernel 

n{t, s) = Xq 9{t — s) / dw/(cj; X , X ) sina;(t — s), (4-2) 
Jo 

and the operator-valued stochastic force ^(t) = Xq^o{Xg,t) in the QBM language. 

In the integrals over the frequency above and below we have assumed a high frequency 
finite cutoff A for the quantum field which regularizes the quantum field's retarded Green's 
function from which we obtain the effective equations of motion of the detector [22]. 

The noise kernel ^{t, s) quantifies the two-time correlation of the Langevin forcing term 

iyit,s) = mt),Cis)}) = Xl da;/(w;x x )cosa;(t-s). (4.3) 

Jo 

[From now on we simply denote I{uj) — I{uj; x^, x^).] Here the average is taken with respect 
to the initial state density matrix of the field. 

By introducing the damping kernel 7(i, s) defined by 

^liT, s) = Mg—j{T, s) = 7(t, s), (4.4) 

we can bring (2.6) to the final form 

MqQ{t) + Mq^lQit) + f dTj{t,T)-(^iT)+2Mqj{t,ti)-Q{U)=i{t), (4.5) 
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where Q'^ = + 7^ is the renormahzed frequency. 
The general solution to (4.1) is given by 



Qit) = Qoit) + f ds G{t,s)i{s), (4.6) 

where G{t,s) is the retarded Green's function for (4.1) satisfying 

MqG{t,s)+Mqn'^G{t,s)- [ dTn{t,T)-G{T,s) = 6{t-s). (4.7) 

For here and below all the Fourier components are defined for positive frequency only, 
which corresponds to Fourier transformation of functions defined on t > 0. Accordingly, 
the Fourier space representation of the above Green's function can be written as: 

G{uj) = {Mq {-uj^ + n^) - f-i{uj)y^ (4.8) 
= {-oj^ - iuj-f{u) + n^y^ M-\ (4.9) 

In the late-time stationary limit, because ft{t,s) leads to dissipation of the detector's free 
motion Qo{t), we see that 

(0(w) ^ G{u;)i{uj). (4.10) 
So the elements of the covariance matrix at late times can be written as 

Vqq = m),Q{t)) = r dujG*{uj)-I{uj)-G{u:), (4.11) 
V^P = {P{t),P{t)) = duj Lo'^M^G*{uj) ■ I{uj) -Giio). (4.12) 



t—^00 



and Vq^ vanishes, as can be inferred from VQp{t) = MVQqit) /2 while VQQ{t) approaches 
an asymptotic constant value. By applying the fluctuation-dissipation theorem [28-32], 
one can eliminate the noise kernel in the covariance matrix elements, giving 



^QQ = ^ r du;lm[Giu;)], (4.13) 
^ Jo 

1 

V^P = - duj uj'^Mllm.[G{u)]. (4.14) 
^ Jo 

For detailed derivation leading to the above results, please refer to Appendix A. 

4.2 Entanglement between Detector and Field in Free Space 

Before we compute the late-time detector- field entanglement with a mirror present, let us 
review the late-time detector-field entanglement for a detector in free space for comparison. 

The first step is to regularize the retarded Green's function of the field at the trajectory 
of the detector. For a detector interacting with a massless scalar field in free space, the 
entire system is governed by the same action as Eq. (2.1) but without the 0:3 > restriction. 
The Dirichlet boundary conditions on the field in the z=0 plane are not needed here. 
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Following [22] we obtain 



(4.15) 



where Qr is the renormalized natural frequency which depends on cutoff A and 7^ = Ay/Svr. 
For free space the retarded Green's function for (4.1) is given by 



(4.16) 



where (ir = ~ 7^ from which the late-time covariances can be computed. 
In the weak coupling limit, up to the first order of 7 , one has 



TrOO 

^QQ,free 



2Mqnr 



TAOO 

^PP,hee 



Mr. 



2 



TT 



:7« 



2(lnA-lnaO 



A2 



1 



(4.17) 



(4.18) 



and the late-time detector-field entanglement in free space is, to the same order of pertur- 
bation in 7q, 



free 



2(lnA-lnf]r) -^-2 
A^ 



(4.19) 



Notice that there is a logarithmic dependence on cutoff scale A, because physically A sets 
a restriction on the ultraviolet modes that the detector can be entangled with. 

4.3 Entanglement between Detector and Field under Mirror Reflection 

In the presence of a perfect mirror, the entire system is governed by action (2.1). The 
Heisenberg equation of motion for the detector after the same regularization as above is 



{d^ + 2j^dt + nl)Q{t) = -^ 

[Cf.(3.4)] and correspondingly 



eit- L)Q{t-L) + Xg^t,x,), 



M, 



(4.20) 



(4.21) 



Here the last term inside the square brackets shows the difference from the free space 
results. It can be interpreted as the contribution from the detector's mirror image located 
at a vertical distance L/2 behind the mirror. 

For this case the exact late-time covariance matrix becomes 



QQjhalf space 



v°° 

* PP, half space 



1 



dwim 



q JO 



Ma 



TT 



dujlui 



-{uj + i7g)2 + + (27^e^^^/L) 

■ 

-{u + i7„)2 + j]2 + (27 e^'^-f'/L) 



(4.22) 
(4.23) 
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Assuming that the detector is only weakly coupled to the field, we can perturbatively 
expand the above integrals and get 



yOO 
QQ, half space 

TAOO 

•^PPjhalf space 



V n 



'.free + ^^QQ 



+ 0(7' 



(4.24) 
(4.25) 



where the terms ^Vqq and 5V^p represent the first corrections to the covariance matrix 
elements due to the presence of the mirror. Physically keeping only these terms in this per- 
turbative expansion is equivalent to ignoring the multiple refiections between the detector 
and the mirror. The exact form for the leading order correction is given below: 



^Vqq 



1 



ttM, 



q JO 
1 1 



dwlm 



Re 



1 



vr MqVtr L 



1 L 



f]2 



(4.26) 



M 

IT 







UJ 



2 / e'^^/L 



(w + i7,)2 + 1^2 \^ - (a; + i7g)2 + r?2 ^ 
Re [{-i + Lrir)) e'^^^r[0, ia,L]] , (4.27) 



in the limit of large cutoff A. 

The change of linear entropy due to the presence of the mirror, compared to the case 
of free space, is given as 

ASl = 5i, half space " ^L.free ~ [e^^'-^r[0, (4.28) 

up to the leading order. Note that despite of the dependence of both Sl, half space and Sl^ free 
on the ultraviolet cutoff A, the difference between them is insensitive to A. 

4.4 Behavior and Physical Interpretation of Late-time Detector-Field Entan- 
glement 

The leading order corrections ASl given in (4.28) is negative (see figure 2 (upper-left)), 
indicating that the presence of the mirror acts to reduce the linear entropy between the 
detector and the field, thereby causing them to be less entangled. Furthermore, in figure 2 
(upper-right), we see that the numerical result of S/,, half space to all orders is monotonically 
decreasing, though wiggling, as L decreases. Clearly in the late-time stationary limit, the 
detector-field entanglement is suppressed in the presence of the mirror compared to the 
case of free space. Such behavior can be understood as follows. 

Since the combined detector-field system is in a global pure state, the detector-field 
entanglement roughly measures the strength of quantum correlation between the two par- 
ties. Indeed, as the detector becomes more and more strongly coupled with the field, 
the correlation between the detector and the field get stronger, and so the detector-field 
entanglement increases monotonically, as shown in the lower plot of Figure 2. 
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Figure 2: (Upper- left) Leading-order corrections to linear entropy given in (4.28) as a 
function of distance between the detector and its image due to the presence of the mirror. 
(Upper-right) Numerical result of the linear entropy of the detector to all orders as a 
function of distance between the detector and its image according to Eqs (4.22) and (4.23). 
(Lower) Linear entropy as a function of 7^. Here Mq = 1, 7^ = 0.02, 17^ = 5. 

The Heisenberg equation of motion (4.20) for the detector's internal degree of freedom 
shows that, through field propagation in the form of spherical wave reflected by the mirror, 
the detector's internal degree of freedom will get a retarded influence from itself after time 
L. There the negative sign attached to the reflected component comes from the Dirichlet 
boundary condition enforced by the mirror. This condition acts on the detector as a time- 
delayed negative feedback. As a result, the correlation established between the detector 
and the field arising from their interactions will be effectively reduced, causing the detector 
to be less entangled with the field. At a closer separation between the detector and the 
mirror, the magnitude of the reflected negative influence becomes larger, causing even 
stronger suppression of the quantum entanglement between the detector and the field. 

4.5 Effect of Mirror Image versus Effect of Real Object 

It is of interest to ask, in our model, whether the effect of the mirror on late-time detector- 
field entanglement can be reproduced by a real detector located at the position of the mirror 
image. In the setup above, the lowest order correction to self-correlators corresponds to the 
physical process in which the detector emits a quanta and then interact with it after it is 
reflected by the mirror. This process contributes a correction of order 0(7g/L), according 
to (4.26) and (4.27). On the other hand, in the cases of two inertial detectors with coupling 
constant equal in magnitude but opposite in signs, as stated in Appendix B, the lowest 
order correction to self-correlators corresponds to the following physical process: detector 
A emits a quanta, which interacts with detector B, then the back reaction from detector B 
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to the field echoes back and interacts with detector A. The contribution of this process is of 
order 0(7g/L^). Therefore, these two cases correspond to two different physical processes 
involving different orders in the coupling constants. 

5 Discussion 

The goal of this paper is to probe into the effect of boundaries on the quantum entanglement 
between local objects and a quantum field, such as that between an atom and its trap or 
cavity surface. Toward this goal we consider the model of an UD detector interacting 
linearly with a massless scalar quantum field under the boundary condition of a perfect 
mirror. Mathematically, the mirror alters the spectral density of the quantum field by 
imposing a Dirichlet boundary condition on it and subsequently affects the dynamics of 
the detector. We describe the detector- field entanglement, both in its early-time dynamics 
and its late-time stationary behavior. We show that in both cases the behavior can be 
depicted in terms of the infiuence of a mirror image. 

One can see more clearly which parties are being entangled if one considers a mi- 
croscopic model of a mirror such as the one considered by Galley et al [35], where the 
mirror's internal degree of freedom is modeled by an oscillator (called mirosc) with very 
light mass (whereby the quantum fiuctuations of the mirror's internal degree of freedom 
will be suppressed) . In this setup we could then consider three physical degrees of freedom: 
the detector, the field and the mirosc. If one first considers the interaction between the 
mirosc and the field, this would yield in the zero mirosc mass or infinite reflectivity limit 
the modified field configuration derived here. Then from the covariance matrix of the de- 
tector one can derive the entanglement between the detector and the modified field. With 
a microphysical model of the mirror one can calculate how the dynamics of the mirosc is 
altered while interacting with the field, how the field is modified, and how the detector 
is entangled with the modified field. The advantage of this approach is that one can see 
clearly how this entanglement can be interpreted as the entanglement between the detector 
and the mirosc (see, e.g., [36] for an example of the successive levels of coarse-graining). 
This reminds us of a similar procedure in electrostatics, namely, how the force between 
a charge and the induced surface charge density on a conducting plate can be calculated 
using the image charge method. 

Another observation using elementary physics of wave reflection upon a mirror is the 
following: Since the field configuration is at the base of inquires into boundary effects on 
the detector field entanglement, the behavior of reflected waves could provide some useful 
guide in building up our intuition on quantum entanglement in this setting. Instead of a 
mirror with perfect reflectivity we can think of two adjoining dielectric media 1, 2 with 
dielectric coefficients ei < £2 (the mirror situation considered above corresponds to the 
case where ei is the vacuum eo, much smaller than €2)- For waves propagating from a 
soft medium 1 to a hard medium 2, the refiected wave is inverted. This shows up in the 
reflected fleld configuration carrying an opposite sign from the original field configuration, 
thus partially canceling it. This cancellation effect is more severe near the mirror surface 
and hence we see the decrease of entanglement as the detector gets closer to the mirror. 
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If this reasoning is correct then in the reverse situation, say a detector in water facing 
the sky (the air is medium 1, where we have assumed £2 > ei), as waves from the heavy 
medium entering a hght medium wih be reflected at the interface with a positive amphtude, 
the entanglement would increase as the detector comes up close to the water-air surface. 
When the experimental techniques improve to the degree that one can measure the quantum 
entanglement between an atom and the trap surface these results could show their practical 
value. 

A natural question to ask is what will be altered if we replace the perfectly reflecting 
mirror here by a more realistic dielectric medium. This is carried out in a sequel paper 
[37] which treats the quantum entanglement between an atom and a dielectric medium by 
means of the influence functional method. As a small corollary from this work we will 
be able to check on the correctness of the above qualitative argument based on symmetry 
and parity considerations. Later papers in this series will address entanglement domain 
and entanglement pattern, and a parallel series on quantum entanglement in topologically 
non-trivial spaces starting with x [38], which can be applied to atoms in a toroidal 
trap. 
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A Derivation of late-time covariance matrix 

According to our definition, the full-time, exact expression for the QQ-element of late time 
covariance matrix is 

VQQ{t)= / doo-I{uj) / dn / dr2G'(ri)cos[^(ri-r2)]G(r2) (A.l) 
Jo Jo Jo 

roo rt rT2+t 

= duj-I{uj) dT2 dfG{f -T2)cos[uj{f -2t2)]G{t2). (A.2) 

Jo Jo Jt2 

At late times the integral can be approximated by 

fOO ft ft 

VQQit)- / duj-I{uj) / dT2 / (ifG'(r-T2)cos[a;(f -2r2)]G'(r2) (A.3) 

Jo Jo Jt2 

roo rt 



= / du-Iiu) df dT2G{f - T2) cos[uj{f - 2t2)]G{t2) (A.4) 
Jo Jo Jo 

= duj- I{uj) drRe | [e-*"^G'(r)] * [e*"^G'(T)] } . (A.5) 

Here * denotes convolution. 

When t — )• -|-oo, the above approximated expressions become exact, and we have 

V^Q= [ du;G*iu)-Iico)-Gico), (A.6) 
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and similarly, the PP-element of the exact late-time covariance matrix is 

V^P = j doj u?MlG*{uj) ■ I{uj) ■ G{uj). (A.7) 

According to previous definitions (4.2) and (4.4) we have 

I{uj) = -ujMgRe[-f{uj)], (A.8) 
vr 

which follows from the fluctuation-dissipation relation between noise and dissipation ker- 
nels, or equivalently their common relation with the spectral density. Then from (4.9), 

G*{uj)-Re[-f{oo)]-G{uj) 

^ \g*{uj) ■ j*{uj) ■ G{uj) + G*{u}) ■ 7(w) • G{u}) 



2ujMq 



V 2iujMq {-2iu!Mg 
Im[G{uj)]. (A.9) 



Therefore in the late-time limit, we can express the covariance matrix elements as (4.13) 
and (4.14) where the explicit reference to the noise kernel is eliminated as G{u;) in (4.8) 
depends only on ^(w). 

B Comparison with the case of two inertial detectors 

In Section VI of Ref. [22], we have obtained the late-time correlators in the case with two 
identical UD detectors at rest at X3 = ±L/2. Let Aq — )• —Xq for the left detector (Qa) and 
Ao — ?• +Ag for the right detector (Qb), one may wonder whether the detector on the right 
{Qb at X3 = +L/2) in this two-detector case would behave the same as the detector at 
the same position (say, Qb at X3 = +L/2) in the above single-detector case with its image 
detector. 

The answer is no. The presence of the other detector separated in a distance L from 
one detector introduces corrections to the late-time correlators of a single detector, which 
are 0{'^'^/ L?) for the self correlators and 0(7g/L) for the cross correlators. This is different 
from the above self correlators Vqq and Vpp, which have 5Vqq and 5Vpp in 0(7g/L). 

This can be understood as follows. The late-time behavior of the correlators in [22] 
are determined by mode functions q^^\t, k) and q^j^\t, k), whose equation of motion reads 
(Eq.(13) in [22] with d and Aq modified) 

{d! + 2^gdt + nl)q^^\t,k) = -?2i0(f-L)gi+)(t-L,k) + Age— *+*^3L/2^ 

{d! + 2^gdt + nl)q^^\t,k) = -^0(t-L)4+)(t-L,k)-Age-"3L/2_ (g^) 

Similar to Eqs. {A6)-(A11) in [22], these equations give the same c^forgjj^^ like (All), 
while the counterpart for Qa is — c^- So the late-time self correlators are still given by Eqs. 
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(48) and (50) in [22], and the cross correlators are those in Eqs. (49) and (51) multiphed by 
— 1. From Eqs. (54) and (55) in [22], one can see that in weak couphng hmit the late-time 
cross correlators are 0(7g/L) and the correction to the self correlators is 0(7q/L^). 



On the other hand, the detector in this paper has 



L 



0{t - L)q+{t - L,k) + 



*sin^, 
M„ 2 ' 



(B.3) 



where q+ (t 
{B.3) gives 



L, k) on the right hand side can be interpreted as the image of g+. Indeed, 



q+{t,k) 



so at late times, 



t»l/75 



-Age 



-iu)t 



sm 



(27„e-^/L) - 



(B.4) 



(Q(i),QW)v = n 



t—^ca 



(27r)32a; 

h 



-(i,k)|^ 



vrM, 



duj Im 



q JO 



1 



— 2i7 w — ilj. 



(27„e--^/L) 



(B.5) 



which is exactly the Vqq in (4.22) after letting h = 1. 

Note that the first order correction to Vqq in free space has exactly the same value as 
the late-time cross correlator {{Qa, Qb}) in the case with two inertial detectors considered 
above. However, the latter will not enter the reduced density matrix of Qb- Compare 



{B.3) with {B.l) and {B.2), one can see that it is {q^^^ + q^^^)/{2i) rather than q^^' has 
the same late-time behavior as q.^ for Mq = 1. 
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